Unusual magnetoresistance in a topological insulator with a single 

ferromagnetic barrier 

B. D. Kong/ Y. G. Semenov/ C. Krowne,^ and K. W. Kim^ 

^Department of Electrical and Computer Engineering, 
North Carolina State University, Raleigh, North Carolina 27695-7911 
^U.S. Naval Research Laboratory, Washington, DC 20375 

Abstract 

Tunneling surface current through a thin ferromagnetic barrier on a three-dimensional topological 
insulator is shown to possess an extraordinary response to the orientation of barrier magnetization. 
In contrast to conventional magnetoresistance devices that are sensitive to the relative alignment of 
two magnetic layers, a drastic change in the transmission current is achieved by a single layer when 
its magnetization rotates by 90 degrees. Numerical estimations predict a giant magnetoresistance 
as large as 800 % at room temperature with the proximate exchange energy of 40 meV at the barrier 
interface. When coupled with electrical control of magnetization direction, this phenomenon may be 
used to enhance the gating function with potentially sharp turn-on/off for low power applications. 
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A topological insulator (TI) is a material that gapless fermion states on its surface and in- 
version electronic bandgap in the bulk are conditioned by exceptionally strong spin-orbital 
interaction.^'- Recent experimental evidences of such surface states in three dimensional 
materials" have sparked significant interests on the subject from the point of view of both 
fundamental physics and potential benefits in numerous device applications. In the imme- 
diate vicinity of the charge neutrality point, the surface states represent an odd number 
of Dirac cones, in which electrons can be described as massless Dirac fermions. Moreover, 
the contact points of the conduction and valence bands are protected by the time reversal 
symmetry that naturally implies suppression of back scattering by nonmagnetic scatter- 
ers of arbitrary potential shape including the point defects. This attribute fundamentally 
distinguishes the TIs from another material with conic band structure, graphene. 

For practical applications, an efficient way of controlling the surface states is crucial. 
It is shown theoretically that TIs can have many exotic physical properties related to the 
breaking of time reversal symmetry. ^'- Of the possibilities, introduction of ferromagnetic ma- 
terials is considered one particularly promising candidate. The exchange interaction with a 
proximate magnetic film can affect the TI surface states far more effectively than externally 
applied magnetic fields. The state-of-the-art synthesis and fabrication technology available 
for magnetic materials/structures also offers nanoscale controllability. Motivated by this 
understanding, several studies have very recently examined their impact on the electron 
transport properties in TIs. In most cases, multiple magnetic layers (or strips) were envi- 
sioned on a TI surface in a manner analogous to the conventional magnetoresistance devices. 
It is not surprising to learn that the modification of electron energy structures induced by 
a combination of disparate magnetization directions can lead to a nontrivial response in 
electrical conduction. One drawback of this scenario, however, is that the functionality of 
the structure is rather limited to realize operations such as logic switching. 

In this study, we exploit a different approach to the current control in a TI with single 
conic surface states. Taking into account that the spin and momentum variables are locked, 
electron transmission through a single ferromagnetic barrier can be significantly dependent 
on the mutual orientation of the incident current J (or, the applied electric field) and the 
magnetization M (see, for example. Fig. 1). For instance, in-plane magnetization of the 
ferromagnetic layer induces a shift in the Dirac cone of the surface states in contact away 
from the Brillouin zone center. When the displacement is in the direction of J, the intro- 
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duced magnetic barrier is not expected to affect electron transmission substantially as the 
conservation of transverse momentum can be maintained. In the general cases, however, the 
momentum conservation rule becomes restrictive, obstructing electron penetration in the 
barrier area. Moreover, orientating M perpendicular to the TI surface leads to a non-zero 
bandgap, i.e., an additional impediment to electron transmission. Thus, a single ferromag- 
netic barrier with variable magnetization appears to be sufficient to modulate the electrical 
current on the surface of a TI. 

The specific structure under consideration is schematically shown in Fig. 1. The TI 
surface is parallel to the x-y plane and a thin strip of ferromagnetic insulator or dielectric 
(FMI) is deposited on top along the y axis. The FMI strip stretches from x = and x = d, 
a width of d that is shorter than the electron mean free path, while its y dimension is 
assumed to be much larger. The current flow is designed to be in the x direction on the TI 
surface. The envisioned magnetoresistance effect can be analyzed, to sufficient accuracy, in 
terms of the surface state Hamiltonian in the lowest order of momentum operator, which 
corresponds to the conic dispersion relation. The quadratic^" and cubic terms responsible 
for the hexagonal warping effects^^ are not taken into account as the interest is primarily 
on the low energy electrons near the Dirac point (i.e., small E). Accordingly the effective 
Hamiltonian is expressed as H = hvp^crxVy ~ o'yPx) + Vg{x) + Hex{x), where p, vp and ai 
are the in-plane wave vector of the electron, its Fermi velocity, and spin operators. Vg{x) 
accounts for the electrostatic potential that may be applied through a gate electrode built 
on the FMI strip, causing Vg{x) = V ioi < x < d and Vg{x) = elsewhere. In the 
mean-field approximation, the exchange interaction Hamiltonian with the proximate FMI 
takes the form H^xix) = q;(x)M ■ a, where the factor a{x) is proportional to the exchange 
integral in the surface region in contact {0 < x < d) and zero otherwise. Combining Hex{x) 
with the first term results in simple momentum renormalization ky ^ ky + A'^^{x) and 
kx — Ay^{x), where notations k = hvpp and A*^^ = a(x)M are used. In this study, 

the magnetic field induced by the FMI strip is not taken into account since its impact 
(typically ~ 0.1 meV in units of energy) is negligible compared to the exchange interaction. 
Thus, it essentially resembles a one-dimensional potential barrier problem with the solutions 
representing the incident (ipi), barrier {ips) and transmitted {ipr) regions. 

In the case of in-plane directed M, the electron transmission problem (with energy E) is 
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solved in the form: 
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where conservation of energy imposes E"^ = k"^ + ky for tpj and tpT, and {E — V)"^ = {k^ — 
Ay^)"^ + {ky + A^^)"^ for ipB- Matching the ip functions at x = and x = d plays the 
role of boundary conditions, defining coefficients t, r, a, and h. Then, the transmission 
probability Pt = |tp is found for an electron as a function of energy E and incident angle. 
A similar procedure can be applied to the case of Mz 7^ 0. Under the ballistic transport 
assumption, the conductance (per unit length in the y dimension) through the FMI barriers 
at temperature T is calculated from the Landauer-Biittiker formalism, '-' 

Gm,{Ef,V) = ^ r dE r deV^{Ep,e,V)F{E,EF,T), (2) 

where F{E, Ef,T) = —df{E,EF,T)/dE and f{E,EF,T) is the Fermi-Dirac distribution 
with Fermi energy Ef- k^ and ky (thus, k^ in region B) are substituted by \E\ cos 6' and 
|i?|sin6' with incident angle 6. The Dirac point is chosen as the reference (zero) for the 
energy coordinate (i.e., Ef, E) in the the present convention. 

Based on this equation [Eq. (2)], the surface resistance Rmj is found by the reciprocal 
relation, Rmj = G^m,? analyzed for three different magnetizations M, each 

directed along the axes j = x, y and z. Specifically, a 30-nm FMI strip is considered at 
room temperature as a function of Ef and V. The exchange energy at the TI/FMI interface 
is assumed to be 40 meV, which appears practically attainable in the current technology. ^'^ 
The calculation results clearly indicate that RMy (j = y) has the smallest resistance of the 
three for all values of Ef and V under investigation. Apparently the Dirac cone displaced 
in the x direction (induced by My) does not significantly hinder electron transmission in 
the same direction in accord with the earlier qualitative analysis. On the other hand, Rm^ 



(j = z) is the largest as expected (due to the bandgap, etc.). As Ep deviates from the Dirac 
point, the filtering effect of the FMI barrier appears to be reduced, leading to a gradual 
decrease in Rm^ and Rm^, while i?My remains unaffected. This can be understood from 
the fact that the current contribution comes predominantly from electrons near the Fermi 
surface (i.e., where df/dE ^ 0). Consequently, an increase of Ep means a higher average 
energy for conducting electrons which are less affected by the "potential barriers". Thus, 
the spin filtering effect is the largest when Ep is at the Dirac point. As for the dependence 
on the gate bias, it is found that RMy increases with V slightly. On the other hand, Rm^ 
and Rm^ decrease as V becomes larger due to the Klein tunneling as discussed in Refs. 4 
and S. In fact, R^iy should also start to decrease at large V (that is beyond the range 
of the current calculation). The values of -Rm^, while they are close, stay generally higher 
than Rm^ as indicated earlier. Since we consider the ideal case where the conduction and 
valence bands are symmetric, the resistance variations described above would also show a 
symmetric behavior for negative values of Ep and V; they are expected to be asymmetric 
in real materials due to the nonlinearity. 

Since the absolute values of the resistance are varied widely by specific parameters for 
materials and structures, the comparison between different magnetic configurations may 
be best achieved by their relative ratios as illustrated in Fig. 2. This is particularly the 
case when one of the quantities (i.e., RMy) is not very sensitive to both Ep and V in the 
range under consideration. As the Fermi energy level increases from to 50 meV, RM^/RMy 
drops from approximately 9.1 to 3.2 while RM^/Ruy is reduced from 6.9 to 3.0. Similarly, 
Rh'ijRMy decreases from 9.1 to 1.8 and Rm^/Rm^ from 6.9 to 1.3 with a gate bias sweep in 
the same range. Adopting the definition of magnetoresistance MR = {RmJ Rmj ~1) x 100%, 
the results clearly suggest that MR as large as 800 % may be achieved at room temperature 
in the proposed single FMI barrier-TI system with a 90° rotation of magnetization [e.g., 
M\\y ^ M\\z, see Fig. 1(a)]. 

Note that the giant magnetoresistance is achieved with only one FMI barrier. This is pos- 
sible because the electron spin and momentum are locked in the manner of Bychkov-Rashba 
effect with an extremely strong spin-orbital interaction — one of the fascinating physical 
phenomena that a TI can provide. In conventional materials, it is necessary to set a refer- 
ence magnetization for spin-based current manipulation. As such, the electric properties are 
varied by controlling the relative spin states between two or more magnetization directions. 
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Ferromagnetic metallic contacts are typically used for spin polarized carrier injection, and 
the current is controlled by another ferromagnetic metallic contact or gate. In a Tl system, 
on the contrary, the current direction itself can work as the reference for spin, eliminating 
the need for one of the magnetic layers. 

For room temperature operation of the structure, the gating function should be robust 
against the thermal broadening, which could wash out the effects of the FMI barrier. Accord- 
ingly, we examine the magnetoresistance as a function of exchange energy and temperature 
as summarized in Fig. 3. Here, two sets of curves are plotted in terms of the exchange energy 
A^^ for analysis. One set represents MR at fixed temperatures (400, 300, 200, 100 K), while 
the other (thin dotted lines) consider it when both temperature and A'^^ vary but their ratio 
is kept constant (i.e., A^^ /ksT = n). Generally speaking, MR increases with A'^^ while 
the temperature has a negative influence as expected. Interestingly, it appears that MR of 
200 % or larger is attainable by the M||?/-M||2; rotation once A^^ satisfies two conditions 
[Fig. 3(a)]; i.e., A^^ is larger than the average electron thermal energy and, at the same 
time, larger than the absolute value of 20 meV. The magnitude of MR increases to over 
1000 %, if is larger than twice the average thermal energy and larger than 20 meV. 
For instance, MR reaches approximately 1300 % when T = 300 K and A'^^ = 52 meV. As 
for the M|||/-M||x rotation [Fig. 3(b)], the values of MR are generally lower than those of 
M||y-M||2; rotation in accord with the earlier discussion. However, MR of a few hundred 
percent can be obtained similarly when A"^^ > ksT. This provides a clear indication that 
the predicted giant magnetoresistance effect could be sufficiently robust for practical use at 
room temperature. 

The induced change in the channel resistance as a function of the magnetization direction 
can find an immediate application in non-volatile memory. There are FMI materials with 
multiple metastable magnetization states (i.e., directions) that can be used for information 
storage. The surface states of a TI in contact would provide an efficient read-out medium. 
Furthermore, one can envision a transistor-like function if dynamic control of magnetization 
direction is achieved as a part of the gate structure. Such a concept may be realized by 
employing a multiferroic heterostructure;^^"^'' that is, an articulated combination of ferro- 
magnetic and piezo-/ferroelectric materials as the gate dielectric. Then, the applied gate 
bias can modulate not only the electrostatic potential in the channel but also the surface 
state band structure itself (i.e., the E — k relationship), leading potentially to a more drastic 
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response in the drain current than that in conventional field-effect transistors (FETs). The 
desired case is when the 90° rotation of magnetization (M||2;— j-MHy) is designed to coin- 
cide with the electrostatic turn-on of the device as schematically illustrated in Fig. 1(c). 
Although the improvement in the on/off slope could be substantial in the scaled devices 
based on tunneling, the largest effect is expected for the intermediate channel lengths where 
the devices switch between the diffusive (low bias) and the ballistic (high bias) transport 
regimes. Due to the break in the time reversal symmetry, it is very likely that the Nl\\z 
case (with a parabolic band structure) possesses a much shorter electron mean free path 
than the other with a massless Dirac cone. Hence, a judicious choice of the channel length 
could ensure a very significant change in the channel mobility in the on/off transition — 
a highly sought-after feature for low-power devices. Quantitative estimate of the potential 
advantage over the conventional FETs requires a detailed investigation of electron-phonon 
interactions^ ' in TIs and is beyond the scope of the current analysis. 

This work was supported, in part, by the SRC Focus Center on Functional Engineered 
Nano Architectonics (FENA) and the US Army Research Office. 
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Figure Captions 



Figure 1. (Color online) Schematic illustration of the considered TI/FMI system. Rotation 
of the magnetization direction around the x (M||i/ ^ M||2;) and z (M||y ^ M||x) axes is 
shown in (a) and (b), respectively, (c) Envisioned device operation of the proposed structure 
as an FET (drain current vs. gate bias). 

Figure 2. (Color online) Relative ratio of surface resistance as a function of (a) Fermi energy 
and (b) gate bias when the magnetization rotates 90° (M||y o M||{x, 2;}). The solid lines 
represent the values of RujRMy [Fig. 1(a)], while the dashed lines show those of RM^/Ruy 
[Fig. 1(b)]. A 30-nm FMI strip is considered at room temperature along with the exchange 
energy of 40 meV at the TI/FMI interface. 

Figure 3. (Color online) Magnetoresistance MR [= {RmJRmj — 1) x 100%] as a function 
of exchange energy A'^^ and temperature T. (a) and (b) correspond to {i=z,j=y) and 
{i=x,j=y), respectively. For both (a) and (b), thick dashed, sohd, dotted, and dash-dotted 
lines represent the results at constant temperature T = 400, 300, 200, 100 K. On the 
other hand, thin dotted curves are obtained with the A'^^ to T ratio kept constant (i.e., 
A^^ /ksT = n, where n = 1, 2, and 3 as denoted in the figure). 



9 




FIG. 1: Kong et al. 
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FIG. 2: Kong et al. 
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FIG. 3: Kong et al. 
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